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ABSTRACT 

We present the current status of cosmic shear based on all surveys done so far. Taken together, they cover more 
about 70 deg 2 and concern more than 3 million galaxies with accurate shape measurement. Theoretical expec- 
tations, observational results and their cosmological interpretations are discussed in the framework of standard 
cosmology and CDM scenarios. The potentials of the next generation cosmic shear surveys are discussed. 
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1. INTRODUCTION 

Gravitational lensing produces distortion of light beams which modifies the image shapes of background galaxies. 
In a Friedman-Robertson- Walker metric, and for stationary, weak gravitational fields, the deflection angle writes 
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V ± $ dl , (1) 



where c is the celerity and <& the 3-dimension gravitational potential. In general, the deflection angle is small 
and lenses can be approximated as thin gravitational systems. This simplifies the relation between the source 
(S) and image (I) positions according to the simple geometrical "lens equation": 

e 1 = 9 s + ^&(0 J ) , (2) 
Uos 

where Dij are angular diameter distances. 

Equations (1) and (2) express how lens properties depend on dark matter distribution and on cosmological 
models. From an observational point of view, gravitational lenses manifest as image multiplication of galaxies 
or quasars, strong and weak distortions of galaxy shape or transient micro-lensing effects. These effects, as 
well as the time delays attached to image multiplications are exploited in order to probe the geometry and 
matter/energy content of the Universe or to observe high-redshift galaxies (see 1, 2, 3, 4 for reviews) 
In the case of weak gravitational lensing, useful approximations relating observed shapes of galaxies to gravi- 
tational shear can be made. In general, the image magnification is characterized by the convergence, k, and by 
the shear components 71 and 72: 

K = \ + ^22); 71 (0) = \ (<PAi ~ ^,22) ; 7a(0) = ¥>,12 = ¥>,2i (3) 
where the (f^j are the second derivatives of <p with respect to the i,j coordinates and 
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The shear applied to lensed galaxies increases their ellipticity along a direction perpendicular to the gradient of 
the projected potential. The lens-induced distortion 6 can then be evaluated from the shape of galaxies as it 

f / (0) 9- 9 d^9 

can be observed from the components of their surface brightness second moment My = - — ' J — : 
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In the weak lensing regime, the relation between the distortion and the gravitational shear simplifies to d « 2^, 
so in principle ellipticity of galaxies, as measured from the second moment components, provides a direct 
estimate of the gravitational shear. The "shear map" can then be used to reconstruct the "mass map" at any 
galaxy position. However, since each galaxy has its own intrinsic ellipticity, and since the background galaxies 
only sparsely sample the sky, the final shear-induced ellipticity map is contaminated by shot noise and cannot 
have infinite resolution. 

The use of ellipticities of galaxies for weak lensing statistics has three important applications: galaxy-galaxy 
lensing, mass reconstruction of clusters of galaxies and gravitational distortion induced by large-scale structures 
of the universe. Each permits to recover properties of the dark matter located in gravitational systems as well as 
the mass density of the universe. In the following we only focus on weak lensing analysis applied to large-scale 
structures, namely the cosmic shear. 

2. MOTIVATIONS FOR COSMIC SHEAR SURVEYS 

Cosmic shear refers to the weak distortion of light bundles produced by the cumulative effects of mass inho- 
mogeneities in the universe. Although theoretical cosmological weak lensing studies started more than 30 years 
ago, it is only recently that this topic moved to an observational cosmology tool. Progresses in the field are 
going incredibly fast, despite observational and technical difficulties, and put cosmic shear as CMB was only few 
years ago. These remarkable progresses were possible because wide field surveys with panoramic CCD cameras 
grow rapidly and can be handled with present-day computing facilities. 

Light propagation through an inhomogencous universe accumulates weak lensing effects over Gigaparsec dis- 
tances. Assuming structures formed from gravitational growth of Gaussian fluctuations, the shape and ampli- 
tude of cosmological weak lensing as function of angular scale can be predicted from Perturbation Theory. To 
first order, the convergence k(9) at angular position 6 is given by the line-of-sight integral 

3 f Zs f x ^ D(z 7)0(7) 
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where \i z ) is the radial distance out to redshift z, D the angular diameter distances, n{z s ) is the redshift 
distribution of the sources. <5 is the mass density contrast responsible for the deflection at redshift z. Its 
amplitude at a given redshift depends on the properties of the (dark) matter power spectrum and its evolution 
with look-back-time. 

The cumulative weak lensing effects of structures induce a shear field which is primarily related to the power 
spectrum of the projected mass density, P K . Its statistical properties can be recovered by the shear top-hat 
variance (10, 19, 21), 

tf) = ^2j o T P K {k)[Ji{k9 c )]\ (7) 

the aperture mass variance (22, 24) 

ooo roc it, 

< M ap> = ^/ o ^P K (k)[-Mk9 c )f, (8) 
and the shear correlation function (10, 19,21): 

^ )e = hJ dk kP.(k)Jo(k9), (9) 



where J n is the Bcssel function of the first kind. P K (k) is directly related to the 3-dimension power spectrum 
of the dark matter along the line of sight, P^d- 



P K (k) = -n 2 m j P 3D [b^, z) W(z, z s ) dz , (10) 

where W(z, z s ) is an efficiency function which depends on the rcdshift distribution of sources and lenses. There- 
fore, in principle an inversion permit to reconstruct the 3-dimension power spectrum of the dark matter from 
the weak distortion field. 

Higher order statistics, like the skewness of the convergence, s 3 (6), can also be computed. They probe non 
Gaussian features in the projected mass density field, like massive clusters or compact groups of galaxies (5' 
6). The amplitude of cosmic shear signal and its sensitivity to cosmology can be illustrated in the fiducial case 
of a power law mass power spectrum with no cosmological constant and a background population at a single 
redshift z. In that case < n{9) 2 > and s 3 (9) write: 

<^) 2 > 1 /^< 7 (^ > l/ 2Ril%(T8fi 0.75 z 0.8^ 2 ) (n) 

and 

S3 W = ^«40n- - 8 ,7 1 - 35 , (12) 

where n is the spectral index of the power spectrum of density fluctuations. Therefore, in principle the de- 
generacy between fi m and as can be broken when both the variance and the skewness of the convergence are 
measured. 



3. DETECTION OF WEAK DISTORTION SIGNAL 
3.1. An observational challenge 

Eq. (11) shows that the amplitude of weak lensing signal is of the order of few percents, which is much 
smaller than the intrinsic dispersion of ellipticity distribution of galaxies. Adopting the null hypothesis that 
no gravitational weak lensing signal is present, we can predict the expected limiting shear if only shot noise 
and sampling are taken into account (in particular, we assume there is no systematic residual errors associated 
to image processing or PSF corrections as those discussed later). For a 3a shear variance limiting detection it 
writes : 

< l(0) 2 >)(Lt= 1-2% 

where At is the total area covered by the survey, a tgal is the intrinsic ellipticity dispersion of galaxy and n 
the galaxy number density of the survey. By comparing with Eq.(ll) one concludes that a cosmic shear survey 
covering 1 deg 2 up to a reasonable depth corresponding to 20 gal./arcmin 2 should in principle be the minimum 
requirements to measure a cosmic shear signal. However, when gravitational shear is included these limits 
change as function of the angular scale and of cosmological parameters. Our simplistic estimate turns out to 
be too pessimistic on small scales because non-linear growth of perturbation producing clusters, groups and 
galaxies amplify the cosmic shear signal by a factor of 2 to 10, depending on cosmologies (5) . Van Waerbeke et 
al (7) explored which strategy would be best suited to probe statistical properties of such a small signal. Their 
results are illustrated on Table 1 and show that the variance of k can already provide a valuable cosmological 
information, even with a survey covering about 1 deg 2 . However, for the skewness one needs at least 10 deg 2 . 
Furthermore, more than 100 deg 2 must be observed in order to uncover information on f^A or the shape of the 
power spectrum over scales larger than 1 degree. 

The shape of objects may also be alterated by atmospheric turbulence, optical and atmospheric distortions 
as well as unexpected telescope oscillations or CCD charge transfer efficiency. Non-cosmological signals reach 
amplitudes one order of magnitude higher than cosmic shear and removing them is indeed the most challenging 
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Table 1. Expected signal-to- noise ratio on the the variance and the skewness of the convergence for two cosmological 
models. In the first column, the size of the field of view (FOV) is given. The signal-to-noise ratio is computed from the 
simulations done by van Waerbeke et al (1999). 
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Table 2. Present status of cosmic shear surveys with published results. 



Telescope 


Pointings 


Total Area 


Lim 


Mag. 


Ref.. 


CFHT 


5 x 30' x30' 


1.7 deg^ 


1= 


=24. 


9[vWME+] 


CTIO 


3 x 40' x40' 


1.5 deg 2 


R= 


=26. 


11[WTK+] 


WHT 


14 x 8' xl5' 


0.5 deg 2 


R= 


=24. 


12[BRE] 


CFHT 


6 x 30' x30' 


1.0 deg 2 


1= 


=24. 


14[KWL] 


VLT/UT1 


50 x 7' x7' 


0.6 deg 2 


1= 


=24. 


15[MvWM+] 


HST/WFPC2 


1 x 4' x42' 


0.05 deg 2 


1= 


=27. 


16 


CFHT 


4 x 120' xl20' 


6.5 deg 2 


1= 


=24. 


17[vWMR+] 


HST/STIS 


121 x 1' xl' 


0.05 deg 2 


V 


* 26 


27 


CFHT 


5 x 126' xl40' 


24. deg 2 


R= 


=23.5 


28 


CFHT 


10 x 126' xl40' 


53. deg 2 


R= 


=23.5 


29 


CFHT 


4 x 120' xl20' 


8.5 deg 2 


1= 


=24. 


30 


HST/WFPC2 


271 x 2.1 x 2.1 


0.36 deg 2 


1= 


23.5 


25 


Keck+WHT 


173 x 2' x 8' 


1.6 deg 2 


R 


=25 


18 




+13 x 16' x 8' 












7 x 16' x 16' 











issues of these surveys. Fortunately, these alterations also modify the shape of stars, hence lensing and non- 
lensing contributions to galaxy shapes can be separated. Over the past decades many works focussed on the 
corrections of non-lcnsing contributions. All follow the basic scheme shown in Figurel, but use different star 
shape analyses and correction techniques (see 3 and 4 for comprehensive reviews). The reliability of artificial 
anisotropy corrections has been discussed and tested at length (8, 9, 12, 13) . In particular, Erben et al (8) 
have demonstrated that the so-called KSB method (23) is able to correct PSF anisotropy very well and allows 
astronomers to detected cosmic shear amplitude down to 1% level with a 10% relative accuracy, even if the 
non-lensing contribution is as large as 10%. 

3.2. Detection 

As we discussed in the previous section, on scale significantly smaller than one degree, non-linear structures 
dominate and increase the amplitude of the lensing signal, making its measurement easier. Few teams started 
such surveys during the past years and succeeded to get a significant signal. Table 2 lists some published results. 
Since each group used different telescopes, adopted different observing strategy and used different data analysis 
techniques, one can figure out their reliability by comparing these surveys. Figure 2 show that they are all in 
very good agreement*. This is a convincing demonstration that the expected correlation of cllipticities is real. 

"The Hoekstra et al (28) data are missing because depth is different so the sources are at lower redshift and the 
amplitude of the shear is not directly comparable to other data plotted 




Figure 1. Basic scheme for the non-lensing distortion correction of cosmic shear surveys. Stars of each field are first 
selected in the vertical branch of the size-magitude plots. The PSF (Raw) is then mapped over the whole field which 
permits to remove the non-lensing distortion anywhere in the field from simple interpolations. The corrected stars provide 
an estimated of the residual (Corrected). 

3.3. Cosmological nature of the signal 

The cosmological origin of the coherent distortion signal detected by all these surveys is not obvious. Even 
if a cosmological signal is present in the data, it could be strongly contaminated by unknown systematic 
contributions or residuals systematic errors produced by the PSF corrections discussed above. An elegant way 
to check whether corrections are correctly done and to confirm the gravitational nature of the signal is to 
decompose the signal into E- and B- modes. The E-mode is a gradient term which contains signal produced 
by gravity-induced distortion. On the other hand, the B-mode is a pure curl-component, so it only contains 
intrinsic ellipticity correlation or systematics residuals. Both modes have been extracted using the aperture 
mass statistics by van Waerbeke et al (17, 26) and Pen et al (30) in the VIRMOS-DESCART" survey as 
well as by Hoekstra et al (29) in the Red Cluster Sequence survey. In both samples, the E-mode dominates 
the signal, although a small residual is detected in the B-mode on small scales for the RCS sample and on 
large scales for the VIRMOS-DESCART sample. These residuals still need further investigations. This strongly 
supports the gravitational origin of the distortion. 

An alternative to a gravitational lensing origin to the signal could be an intrinsic correlations of ellipticities 
of galaxies produced by proximity effects (tidal torques). Several recent numerical and theoretical studies (31, 

f For the VIRMOS spectroscopic survey, see (34). For the VIRMOS-DESCART survey see 

ht t p : / / 1 er apix . iap . fr / D escart / 

*Data of the VIRMOS-DESCART survey were processed and analysed at the TERAPIX data center: 
http://terapix.iap.fr 
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Figure 2. Top hat variance of shear as function of angular scale from 6 cosmic shear surveys. The acronyms of the right 
refer to the references reported in Table 2. The open black stars illustrate the expected performances of the CFHT-LS 
which will start by 2003 with Megacam at CFHT. This is the expected signal from the "Wide Survey" which will cover 
about 200 deg 2 up to Iab = 24.5. As one can see, for most points the errors are smaller than the stars. 

32) have addressed this point. From these first investigations, it seems that intrinsic correlations are negligible 
on scales beyond one arc-minute, provided the survey is deep enough, but this critical point is still debated. 
Nevertheless, we hope that deep surveys are reliable since in that case, most lensed galaxies along a line of sight 
are spread over Gigaparsec scales and have no physical relation with its apparent neighbors. Hence, since most 
cosmic survey are deep, we expect they are almost free of intrinsic correlations. 

4. COSMOLOGICAL INTERPRETATIONS 
4.1. 2-point statistics and variance 

A comparison of the top- hat shear variance with some realistic cosmological models is ploted in Figure 2. Its 
amplitude has been scaled using photo- z which gives < z >~ 1. On this plot, we see that standard COBE- 
normalized CDM is ruled at a 10— a confidence level. However, the degeneracy between Q m and a% discussed 
in the previous section still hampers a strong discrimination among most popular cosmological models. The 
present-day constraints resulting from independent analyses by Maoli et al (15), Rhodes et al (16), van Waerbeke 
et al (17,26), Hoekstra et al (28,29) and Refregier et al (25) can be summarized by the following conservative 
boundaries (90% confidence level): 

0.05 <n m < 0.8 and 0.5 < <r$ < 1.2 , (14) 
and, in the case of a flat-universe with fi m = 0.3, they lead to a$ ~ 0.9. 
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Figure 3. Constraints on Q and ag for the fiat cosmologies. The confidence levels are [68, 95, 99.9] from the brightest 
to the darkest area. The gray area and the dashed contours correspond to the computations with a full marginalisation 
over the default prior F £ [0.05,0.7] and z s £ [0.24,0.64]. The thick solid line contours are obtained from the prior 
T € [O.f, 0.4] and z s € [0.39, 0.54] (which is a mean redshift z s £ [0.8, 1.1]). From van Waerbeke et al. (2002) (26). 

4.2. Toward breaking the Q m — as degeneracy 

The measurement of non-Gaussian features needs informations on higher order statistics than variance. Eq.(12) 
shows that the skewness of k seems the best suited quantity for this purpose. However, in contrast with the 
variance, the skewness demands first to reconstruct the mass map of the field, which complicates the process. In 
fact, its measurements suffers from a number a practical difficulties which are not yet fixed. In particular, the 
masking process generates a Swiss-Cheese pattern over the field which significantly degrades the quality of the 
mass reconstruction and put strong doubts on the reliability of its skewness and its cosmological interpretation. 
Recently, Bernardeau, van Waerbeke & Mellier (37) have proposed an alternative method using some specific 
patterns in the shear three-point function. Despite the complicated shape of the three-point correlation pattern, 
they uncovered it can be used for the measurement of non-Gaussian features. Their detection strategy based on 
their method has been tested on ray tracing simulations and turns out to be robust, usable in patchy catalogs, 
and quite insensitive to the topology of the survey. 

Bernardeau et al (38) used the analysis of the 3-point correlations function on the VIRMOS-DESCART data. 
Their results on Figure 4 show a 2.4a signal over four independent angular bins, or equivalently, a 4.9-cr 
confidence level detection with respect to measurements errors on scale of about 2 to 4 arc-minutes. The 
amplitude and the shape of the signal are consistent with theoretical expectations obtained from ray-tracing 
simulations. This result supports the idea that the measure corresponds to a cosmological signal due to the 
gravitational instability dynamics. Moreover, its properties could be used to put constraints on the cosmological 
parameters, in particular on the density parameter of the Universe. Although the errors are still large to permit 
secure conclusions, one clearly see that the amplitude and the shape of the 3-point correlations function match 
the most likely cosmological models. Remarkably, the ACDM scenario perfectly fit the data points. 
The Bernardeau et al (38) result is the first detection of non-Gaussian features in a cosmic shear survey and 
it opens the route to break the O m — as degeneracy. Furthermore, this method is weakly dependent on other 
parameters, like the cosmological constant or the properties of the power spectrum. However, there are still 
some caveats which may be considered seriously One difficulty is the source clustering which could significantly 
perturb high-order statistics (Hamana et al 2000 (33)). If so, multi-lens plane cosmic shear analysis will be 
necessary which implies a good knowledge of the redshift distribution. For very deep cosmic shear surveys 
which probe galaxies fainter than the limiting spectroscopic capabilities of 10-meter class telescopes, the lack 




Figure 4. Results for the VIRMOS-DESCART survey of the reduced three point correlation function (38). The solid 
line with error bars shows the raw results, when both the E and B contributions to the two-point correlation functions 
are included. The dot-dashed line with error bars corresponds to measurements where the contribution of the B mode 
has been subtracted out from the two-point correlation function. These measurements are compared to results obtained 
in rCDM, OCDM and ACDM simulations (dashed, dotted and dot-dashed lines respectively). 



of reliable photometric redshift calibration be could be a serious concern. 

5. OUTLOOK 

Because on going surveys increase both in solid angle and in number of galaxies, they will quickly improve 
the accuracy of cosmic shear measurements, at a level where fi m and erg will be known with a 10% accuracy. 
Since it is based on gravitational deflection by intervening matter spread over cosmological scales, the shape 
of the distortion field also probes directly the shape of the projected dark matter power spectrum. Pen et 
al (30) explored its properties by measuring for the first time the C(l) of the dark matter (see Figure 5). 
Although the reliability of these C{1) is still unknown, this work shows this is feasible with present data. The 
capabilities of cosmic shear surveys to probe dark matter can be also used to explore the relations between 
light and mass at different angular scale, different redshifts and in different environment. Recently Hoekstra 
et al (39) joined together the RCS and the VIRMOS-DESCART surveys and computed both the biasing and 
the galaxy-mass cross-correlation. They used the aperture mass and aperture number density statistics, as 
it was originally proposed by Schneider (35) and van Waerbeke (36). They found that the biasing is linear 
on scale beyond 3 h , but turns out to be more complex on lower scale, indicating a significant non-linear 
and/or stochastic biasing. These results are indeed difficult to interpret in a cosmological context without high 
resolution numerical simulations and a detailed scenario of galaxy formation in hands. But they show the 
remarkable potential of cosmic shear survey for challenging theory of galaxy formation against observations. 
Moreover, it is worth noticing that this work was made possible only because the two surveys merged their data 
in order to cover a sufficiently large field of view. The investigation of mass and light relations is certainly one 
of the most promising goals of wide cosmic shear surveys in the future. 

However, we expect much more within the next decade. Several new wide field surveys are in progress 
at ESO/La Silla (Erben et al 2002 in preparation), at CTIO (Jarvis et al), at NO AO (Januzzi et al) and at 
SUBARU (Miyazaki 2002 (41)). Surveys covering hundreds of degrees, with multi-bands data in order to get 
redshift of sources and possibly detailed information of their clustering properties, are scheduled. The CFHT 
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Figure 5. Cosmological results from cosmic shear surveys: The angular power spectrum of the convergence field from 
the VIRMOS-DESCART survey is plotted (From Pen et al 2001). These are the first C(l) of dark matter ever measured 
in cosmology. 

Legacy Survey^ will cover 200 deg 2 and is one of those next-generation cosmic shear survey. Figures 1 and 
6 shows their expected potential for cosmology. On figure 1 we simulated the expected signal to noise of the 
shear variance as function of angular scale for a ACDM cosmology. The error bars are considerably reduced 
as compared to present-day survey. On Fig. 6, we compare the expected signal to noise ratio of the CFHT 
Legacy Survey with the expected amplitude of the angular power spectrum for several theoretical quintessence 
fields models. It shows that even with 200 deg 2 which include multi-color informations in order to get redshift 
of sources, one can already obtain interesting constraints on cosmology beyond standard models. 
The use of cosmic shear data can be much more efficient if they are used together with other surveys, like 
CMB (Boomerang, MAP, Planck), SNIa surveys, or even galaxy surveys (2dF, SDSS). For example, SDSS 
will soon provide the 100, 000 quasars with redshifts. Menard & Bartelmann (42) have recently explored the 
interest of this survey in order to cross-correlate the foreground galaxy distribution with the quasar population. 
The expected magnification bias generated by dark matter associated with foreground structures as mapped 
by galaxies depends on fl m and the biasing erg. In principle magnification bias in the SDSS quasar sample 
can provide similar constrains as cosmic shear. Both cosmic shear and cosmic magnification wide field surveys 
should soon produce impressive constraints on cosmological models. 
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